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Abstract 

The cross sections describing the transverse momentum distributions of Higgs bosons 
^ , are discussed within the framework of unintegrated parton distributions obtained from 

QJ \ the CCFM equations in the single loop approximation. It is shown how the approxi- 

' mate treatment of the CCFM equations generates the standard expressions describing 

the soft gluon resummation effects in the corresponding cross sections. Possible differ- 
^ ■ ences between exact and approximate solutions of the CCFM equations are discussed on 

' the example of the gg ^ H fusion mechanism, which gives the dominant contribution to 

Higgs production. 



1 Introduction 



The transverse momentum distributions of Drell-Yan pairs, the W or Zq gauge bosons or the 
Higgs bosons produced in high energy hadronic coUisions are known to be strongly affected by 
the soft gluon (recoil) resummation effects - ^7]. These effects are particular important in 
the region p% Q'^-, where Q and Pt denote the mass and transverse momentum of the virtual 
photon 7* or the Zq or Higgs bosons, respectively. The soft gluon resummation effects are 
also important in the prompt photon production in hadronic collisions. 

All these reactions reflect the basic partonic fusion subprocesses, i.e. g + g — >■ 7* ^ 
q + q ^ W, q + q ^ Zq, g + g ^ H, g + q{q) — > q{q) + 7 etc. The transverse momenta of the 
lepton pairs, the W, Zq and/or Higgs bosons etc. should therefore just be equal to the sum of 
the transverse momenta of the colliding partons. The transverse momentum distributions of 
the lepton pairs, the W, Zq and/or Higgs bosons etc. produced through the fusion of partons 
do therefore probe the transverse momentum distributions of the partons, i.e. the unintegrated 
parton distributions [SI ll],[IH]-[2II. The corresponding cross sections are determined by the 

factorisation - like convolutions (subject to the transverse momentum conservation) of the 
unintegrated parton distributions with the corresponding partonic cross sections. In this simple 
extension of the collinear picture to the more exclusive configuration, in which the transverse 
momenta are not integrated over the gluon resummation effects, are naturally attributed to 
the unintegrated parton distributions. Evolution of the unintegrated distributions is described 
in perturbative QCD by the Catani, Ciafaloni, Fiorani, Marchesini (CCFM) equation j2Hl- It 
corresponds to the sum of the ladder diagrams with the angular ordering along the chain. The 
gluon radiation controlling evolution of the unintegrated parton distributions generates the 
transverse momentum of the partons through the simple recoil effect. The CCFM equation 
interpolates between the (LO) DGLAP evolution at large and moderately small values of x and 
the BFKL dynamics at low x. Pure DGLAP evolution corresponds to the so called "single loop 
approximation" of the CCFM evolution I2H1-E2]- 

The purpose of this paper is to examine the transverse momentum distributions of the Higgs 
bosons produced in hadronic collisions using the unintegrated parton distributions framework. 
To be precise we shall limit ourselves to the gg ^ H fusion mechanism and determine the 
corresponding cross sections in terms of the unintegrated gluon distributions. The latter will 
be obtained from the solution of the CCFM equation in the single loop approximation, which 
may be adequate for the moderately small values of x probed in the (central) Higgs boson 
production. We shall utilise the fact that this equation can be diagonalised by the Fourier- 
Bessel transform and becomes evolution equation for the (scale dependent) parton distributions 
fi{x,b,Q) for fixed b, where 6 is a vector of impact parameter conjugate to the transverse 
momentum. We express the cross sections in terms of the function /G(a^i,25 b, Q) and show that 
the CCFM equation correctly reproduces the double and single ln(Q^6^) resummation in the 
limit 3> 1, that corresponds to <^Q'^ . 

The content of our paper is as follows. In the next section we recall the basic formulas 
describing the soft gluon resummation effects in the transverse momentum distributions de- 
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scribing of Higgs bosons produced through the gluon-gluon fusion mechanism. In section 3 we 
discuss this distribution within the unintegrated gluon distributions framework. We point out 
that the basic quantities controlhng the transverse momentum distributions in this framework 
are the unintegrated distributions in the 6— representation. In section 4 we discuss the CCFM 
equation in this representation. Section 5 contains discussion of the soft gluon resummation 
formulas as the result of the approximate treatment of the solution(s) of the CCFM equation 
in the 6— representation. In section 6 we present numerical results concerning the transverse 
momentum distributions of Higgs bosons and compare predictions based upon exact solution 
of the CCFM equation with the standard LO soft gluon resummation formulas. Finally, in 
section 7 we give summary of our results. 



2 Soft gluon resummation formulas 

The soft gluon resummation effects are most conveniently formulated using the 6-representation. 
Thus the cross-section for Higgs boson production in pp collisions corresponding to the gluon - 
gluon fusion subprocess gg H is expressed in the following way jTH [T7] 



da (To Q 



dQ'^dydp\ rs 



1 

7r2(5(g2-m|)- / dbbJoibpT)WGGixi,X2,b,Q), (1) 
^ Jo 



with Jo{z) being the Bessel function. In eq. pt and y denote the transverse momentum 
and rapidity of the Higgs boson, respectively, and mn is the Higgs boson mass while xi^2 are 
the longitudinal momentum fractions of the gluons given by 



xi = \ exp{y) 



'Q2 


+ pI 


s 






s 



X2 = \l / ^M-y) (2) 

where s is the square of the CM energy of the colliding pp system, and 

r = X1X2 . (3) 

The cross section ctq is given by: 

^0 = ^5^^ asif^l) (4) 
57o TT 

where Gp is the Fermi constant. The corresponding cross-section for the Drell-Yan production 
of muon pairs reads: 



da a^^ 1 



dQ'^dydp\. r 2 



dbbUbpT)Y.e1{W,^,, {Xi,X2, b,Q) + Wq^qXxi,X2,b,Q)} (5) 
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where now p^, Q and y denote the transverse momentum, invariant mass and rapidity of the 
lepton pair, respectively, and 0:1^2 denote the momentum fractions of the parent hadrons carried 
by a quark or antiquark and are defined by eq. (j2I). The elementary DY cross section a^^ is 
given by: 

=W 

The basic functions driving the corresponding cross sections are the functions Wij{xi, X2, b, Q) 
which can be interpreted as the scale dependent parton luminosities in the 6— representation. 
The soft gluon resummation gives the following representation of these functions 

W,,ixi,X2,b,Q) = W^y^(xi,X2,&,Q) exp|i[S,(6,g) + S',(6,g)] 

X x,pt"ix^,cyb^)x2pf^ix2,cyb^). (7) 

The functions p'^^^ are given by the conventional integrated parton distributions Pi probed at 
the scale /z^ = with c ~ 1 0, i.e. 

pt"ix,fi^) = Pi{x,fi') + 0(a,(/i2)) (8) 
with the 0{as{fi'^)) corrections given in [7, 8j. The exponents Si are given by: 



S^{b,Q) 



2/f,2 



Q 



(9) 



where the functions Ai and Bi are defined by the perturbative expansion: 

n=l 

00 

■a. 



BMs) = Bl-'- (10) 

n=l 

The LO coefficients are given by [H]: 

Ag) = 2Ca 

= iNjTn - (11) 

= bP = (12) 

where for the SU{3)c we have: Ca = 3,Cf = 4/3, Tr = 1/2, and Nf is the number of active 
flavours. Higher order coefficients are defined in [TU]. The factors W^^^{xi,X2,b,Q) describe 
the non-perturbative contribution(s) [131 CH HZI ■ 
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Coefficients A\"'^ and -B^-"^ are directly connected to the part of the sphtting functions Pij{z) 
corresponding to real emission. To this aim we introduce the functions Rg{z) and Rq{z) defined 
by: 

Rg{z) = z[P,,iz) + Pq,iz)] (13) 
R,{z) = z[P,,{z) + P,,{z)] , (14) 

and rearrange these functions as below 

Rg{z) = + Rg{z) (15) 

where = Cp and = 2Ca, and Rg{z) and Rq{z) are regular at 2; = 1. In LO we get the 
following relations 

A^^ = r,, i = q,G (17) 

= 2 [ dzRiiz). (18) 
Jo 

It is also useful to remind the following identity: 

/ dz[R,{z)-Pgg{z)] = 0. (19) 
Jo 

3 Transverse momentum distributions of Higgs bosons 
and DY pairs 

The cross section for Higgs production in pp collisions corresponding to the gluon-gluon fusion 
subprocess gg ^ H can be expressed in the following way in terms of the scale dependent 
unintegrated gluon distributions /G(a;i,2, ^t, Q) 



dQ'^dydpj. ts 

I J2l, ^2 

X 



J ik^^^^^^^^^^^^^^^^^^^^^^^^(2)(ki + k2-pT). (20) 

In this equation px and y denote, as before, the transverse momentum and rapidity of the 
Higgs boson, respectively, rriH is the Higgs boson mass, ki^2 are the transverse momenta of 
the colliding gluons and xi,2 are their longitudinal momentum fractions (see eq. (j21))- The 
cross-section (Tq is given by eq. (jH). From eqs. (j2l3|) we also have 

TS = Q'^+pI (21) 
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Substituting in eq. (pUj) the following representation of the 6 function, 



^(2) (k^ + _ p^) = _L_ j d% exp{-z(ki + k2 - PT)b} , (22) 

we find that the integral defining the differential cross section for Higgs production can be 
directly expressed in terms of the unintegrated gluon distributions fG{x,b,Q) in the impact 
parameter representation 

poo 

Ti^SiQ^ - ml) / db^ fcixi, 6, Q) /g(x2, b, Q) Mbpr) . (23) 

^0 



da aoQ 



dQ^dydp^ ts 



where 



POO 

fGix,b,Q) = / dkTkTMbkT)fGix,kT,b). (24) 

^0 



The cross section for Higgs production which is driven by the gg ^ H fusion is determined 
by the gluon distributions. Similar representation of the cross section, but in terms of the 
quark and antiquark distributions in the b space, can be obtained in the case of the Drell-Yan 
production of the muon pairs. The latter reaction reflects the process gg ^ 7* ^ f^~^f^~, and 
the corresponding expression for the cross-section reads 

da <^ 



dQ2QyQp2^ r J db'MbpT)J2^^ [fl{xi,b,Q)fl{x„b,Q) + fl{x„b,Q)fl{x„b,Q)] 

(25) 

where now p^, Q and y denote the transverse momentum, invariant mass and rapidity of the 
lepton pair, respectively, Xi^2 denote the momentum fractions (j21) of the parent hadrons carried 
by a quark or antiquark. The functions f^f and f^f are the quark and antiquark unintegrated 
distributions in the 6— representation. As in the case for the gluons, they are linked to the scale 
dependent fc^ distributions of quarks and antiquarks through the Fourier-Bessel transform (cf. 
eq. (j21|))- The elementary DY cross section (T,f^ is given by eq. ©. Comparing eqs. 

with eqs. (^E)), we find that within the unintegrated parton distrtibution framework the 
functions Wij{xi,X2,b,Q) are given by: 

Wi,ixi,X2,b,Q) = 4/,(xi,6,g)/,(x2,6,g). (26) 



4 Single loop CCFM equation in the 6— representation 

In this section we recall the CCFM equation in the single loop approximation. Then we 
extend this equation to the system of equations incorporating unintegrated quark and antiquark 
distributions, fq. and /g.. The original Catani-Ciafaloni-Fiorani-Marchesini (CCFM) equation 
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|2Hll2ni for the unintegrated, scale-dependent gluon distribution fG{x,kT,Q) has the following 
form: 

fG{x,kT,Q) = faixM.Q) + [' -Q{Q-qz) e{q- go) ^MQ,q,z) 

J Jx z 2tt 

/(|,|kT+(l-z)gU), (27) 

where As{Q, q, z) and /S.Ns{kT, q, z) are the Sudakov and non-Sudakov form factors, 

^s{Q.q.z) = exp{- 4^ / dzzP,,{z) \ (28) 

[ J{qz)^ P Jo J 

AMs{kT,q,z) = exp\- — "V^f^ \ ■ 29 

The variables x, kx and Q denote the longitudinal momentum fraction, the transverse momen- 
tum of the gluon and the hard scale, respectively. Eq. is generated by the sum of ladder 
diagrams with angular ordering along the chain. The hard scale Q is defined in terms of the 
maximum emission angle [211 12H]- The constraint Q{Q — qz) in eq. (jTTj) reflects the angular 
ordering, and the inhomogeneous term, f^{x,kT,Q), is related to the input non-perturbative 
gluon distribution. It also contains effects of both the Sudakov and non-Sudakov form-factors 

m- 

Eq. ()27|) interpolates between a part of the DGLAP evolution at large and moderately small 
values of x, and the BFKL dynamics at small x. To be precise, it contains only the g ^ gg 
splittings and only those parts of the splitting function Pgg{z) which are singular at either 
z — > 1 or z — i> 0. Following |33J|SS] we shall make the following extension of the original CCFM 
equation (jTrj) : 

1. We introduce, besides the unintegrated gluon distributions faix, kx, Q), also the uninte- 
grated quark and antiquark gluon distributions, fqX^, kr, Q) and /^^(x, fcy, Q)- 

2. We include, in addition to the g gg splittings, also the q gq, q gq, and g qq 
transitions along the chain. 

3. We take into account the complete splitting functions Pab{z), and not only their singular 
parts. 

These extensions make the CCFM framework more accurate in the region of large and moder- 
ately small values of x (i.e. x > 0.01 or so). In this region one can also introduce the single-loop 
approximation that corresponds to the following replacements: 

e{Q-qz) e(Q-g) 
ANsikT,q,z) 1. (30) 



X 



2NcANs{kT,q,z) + 



2NrZ 
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It is also convenient to "unfold" the Sudakov form factor such that the real emission and 
virtual terms appear on an equal footing in the corresponding evolution equations. Finally, 
the unfolded system of the extended CCFM equations in the single-loop approximation has the 
following form: 



X Pqq{z) Q{Z-X) fi^s k'T, - f^gi^, ^T, Q) 



(31) 



fs{x,kT,Q) = fs(x,kT) + f'dz f ^,^^e(q'-q'o)m-Q) 

Jo J T^T 27r 



- Pqq{z) fs{x,kT,q) 



Jo J T^T ^TT 

X {q{z-x) [p,,{z)fs[-^,k'^,q)+Pg,{z)fG[-^,k'.,,q) 



- [^Pggiz) + ^-^99(^)1 faix, kr, q) 



(32) 



(33) 



where 

Ic't = kT + (1 - ^)q. (34) 

The functions f^si^j ^t, Q) are the unintegrated non-singlet quark distributions, while the 
unintegrated singlet distribution, fs{x,kT,Q), is defined as 

fs{x, kr, Q) C-^^^^^' ^) + ■ (35) 



1=1 



The functions Pab{z) are the LO splitting functions corresponding to real emissions, i.e.: 

P,,{z) = N, [z' + (1 - zf] 



Pg,{z) = 2iV, 



z 

z 1- z ,^ . 
^ +z{\-z) 

L — Z Z 



(36) 
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where Nj and Nc denote the number of flavours and colours, respectively, and Cp = {N^ — 
l)/{2Nc). It should be observed that after integrating both sides of eqs. (j^ - ljHHj) over cPk-r, 
we get the usual DGLAP equations for the integrated parton distributions pk{x,Q'^), defined 
by 



dk^fkix, kr, Q). 



(37) 



A very important merit of eqs. (j^ - lfH^ is the fact that they can be diagonalised by the 
Fourier- Bessel transform 



(38) 



fk{x,kT,Q) = I dbbJo{bkT)fk{x,b,Q) 
where the inverse relation reads 

fk{x,b,Q) = I dkrkx Joibkr) fk{x,kT,Q) . 



(39) 



In the above k = NS, S, G, and Jo is the Bessel function. At 6 = the functions /fc(x, b, Q) are 
related to the integrated distributions J9j(x, Q^) 



fk{x,b = 0,Q) = -xpk{x,Q^). 



(40) 



The corresponding evolution equations for f^s, fs and fc, which follow from eqs. (p?T|l - (p!?jl . 
have the following form: 



Q'- 



dfNs{x,b,Q) asiQ'^ 



2n 



dzPqq{z) 



Qiz -x) Mil -z)Qb)f MS {^,b,Q) 



zPqq{z) fNs{x,b,Q) 



(41) 



dfsix,b,Q) as{Q 



2\ /•! 



271 



Q{z - x) Mil - z)Qb) 



Pgqiz)fs{-,b,Q 



X 



+ Pqgiz)fG{-,b,Q 



X 



[zP,g{z)+zPg,{z)]fs{x,b,Q) 



dfG{x,b,Q) a,(Q2 



2tt 



e{z - x) Mi'^ - z)Qb) 



PUz)fs{-b,Q 



- [zPggiz) + ZPggiz)] /g(X, b, Q) 



with the initial conditions 



/fc(x,6,Qo) = FkixM 



(42) 



(43) 
(44) 
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We have therefore found that the CCFM scheme in the single loop approximation does directly 
provide the system of evolution equations for the unintegrated parton distributions fk{x,b,Q) 
in the 6— representations. These functions are the basic quantities in the description of pt 
spectra, eqs. within the unintegrated parton distributions framework. In the next 

section we discuss their solution and show that they correctly incorporate the resummation 
effects. 

5 Soft gluon resummation formulas from the solution of 
the CCFM equations 

In order to understand the structure of the solution of the CCFM equations (PT|) - (PH|l . it is 
convenient to introduce the moment functions fi{n,b,Q) 

Mn,b,Q) = [ dxx''-'Mx,b,Q) (45) 

We shall consider at first the non-singlet, valence quark distributions for which the CCFM 
equation in the moment space reads: 

Its solution reads: 

fNs{n,b,Q) = fNs{n,b,Qo) exp{I{n,b,Q)} (47) 
where the function I{n, b, Q) is given by: 

I{n,b,Q)= / 4^7^/ dzP,,{z){z-U{l-z)Qb)-l} , (48) 

JQI T ^TT Jo 

and the function fNs{n,b,Qo) is given by the input unintegrated distribution at the reference 
scale Qo- At 6 = we have: 

fNsin, 6 = 0, Q) = ^ pr,s{n, Q^) (49) 

fNs{n,b = 0,Qo) = ^PNs{n,Ql) (50) 

where p{n, Q"^) are the moment functions of the integrated distributions. We shall now rearrange 
the integral I{n, b, Q) as below: 

/(n,6,Q) = /i(n,6,Q) + h{b,Q) (51) 



10 



where 



h{n,b,Q) = [ ^^All f dz P,,{z) {z- - 1) Mil - z)Qb) (52) 

h{b,Q) = r [\z P,,{z) {Mil - z)Qb) - 1} . (53) 

Jq2 g2 27C Jo 



The dominant contribution to the integrals Ji and I2 comes from the region in which argu- 
ment of the Bessel function Jo{u) is small, where the Bessel function itself can be approximated 
by unity, i.e. we can adopt the following approximation: 

Mu) ~ e(c-M) (54) 

with c being of the order of 1; for simplicity we set c = 1. The dominant contribution to the 
integral Ji comes from the region Qb < 1 and adopting approximation (jMj) we get for b > 1/Q: 

hin,b,Q) ^ f" ^^fifl f dzP,,{z){z--l). (55) 
JQI r 27r Jo 

We observe that approximate representation of /i(n, 6, Q) is equal to the argument of the 
exponent describing evolution of the moments of the integrated non-singlet quark distributions 
from the scale Qq to the scale 1/6^. For 6 -C 1/Qo we may set 

fNsin,b,Qo) ~ fNs{n,b = 0,Qo) = ^PNs{n,Ql) (56) 

to finally get 

fNsin,b,Qo) exp{Ii{n,b,Q)} ~ • (57) 

After inverting moments one finds 

fNs{x,b,Q) ~ ^xpNs{x,l/b'^) . (58) 

Comparing p6|) with ((Zj) we find that in the functions Wij{xi,X2,b,Q) defined by the non- 
singlet quark distributions, we have identified the first factor in the soft gluon resummation 
formula to be equal to the product of the integrated parton distributions at the scale of the 
order of 1/6^. 

We now show that after using approximation in the evaluation of l2{b,Q), the form 
factor exp{2/2(6, Q)} becomes equal to the form factor exp{Sg{b, Q)} in the LO approximation 
with the exponent Sq{b,Q) defined by eq. To this aim we note that after rearranging the 
splitting function Pqq{z) 



^-{1 + z) 



(59) 
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we can represent the integral l2{b, Q) as the sum of two terms: 



h{h,Q) = ll{h,Q) 



• Q) 



(60) 



where 



JqI 27r Jo 



dz 
1 - z 



[J,{{l-z)qh)-l] 



Jq2 g2 27r Jo 



dz{l + z) [Jo{{l - z)qh) - I] 



(61) 
(62) 



Using approximation (jH^ we get: 



ll{h,Q) -Q{Q^ - l/h-)C, I In(gV) 



(63) 



e{Q'-l/b^)-CF 



a,(g2) 



(64) 



Thus ffHUI) becomes 



2/2(6,g) ^ -2e(Q2_i/62)c'p 



0' .7^2 



/62 g2 27r 



ln(g 6 



(65) 



We note that this expression coincides with that given by eq. Q defing the exponent Sq{h, Q) in 
the LO approximation. To be precise, the first term corresponding to I^ib, Q) is identically equal 
to the term defined by Aq{as) for the fixed coupling only and there are subleading differences 
due to different scale of the QCD coupling. It may be observed, however, that we would get 
exactly the same result for this term choosing the scale of the coupling equal to (g(l — z)Y 
instead of q^ in the integrals defining l\{h,Q). The second term corresponding to /|(6, Q) is 
identical to that defined by Bq{as)- Finally we observe that the non-perturbative contribution 
factorises as the Q independent factor, see eq. (jUj). 

The valence quarks dominate in the production of DY muon pairs and in the production 
of gauge or bosons in pp collisions. However, inclusive production of the Higgs bosons, 
which is dominated by the gluon-gluon fusion, requires the knowledge of gluon and sea quark 
distributions. We shall now derive the soft gluon resummation formula for the unintegrated 
gluon and sea quark distributions, and for the function WGG{.xi,X2,h,Q), starting from the 
system of the CCFM equations ()42I43|) . In this case we cannot solve the system of the CCFM 
equations in an analytic form and the soft gluon resummation expressions will follow from the 
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approximate treatment of the CCFM equations themselves. To this aim it is convenient to 
integrate these equations and represent them as the system of integral equations: 



fs{x,b,Q) = Ps{x,b)+ / dz 



g2 27r 



0(g' - qt) 0(Q - q) 



X \ Q{z-x)U{l-z)qh) P,,iz)fs[^,b,q)+P,g{z)fG[^,b,q 



[zPgg{z) + ZPgq{z)] f g {X , 6, q) 



(66) 



/g(x,6,Q) = Pa{x,b)+I\z I ^^^e{q'-ql)e{Q~q) 

X |e(z - x) Mil - z)qb) [Pggiz) fs b, g) + P,,(z) fa b, q 



- [zPgg{z) + zP,g{z)] faix, b,q)y (67) 
We rearrange these equations by adding and subtracting in the integrand of (jMj) the term 



27rg2 



[zPgg{z) + zPgg{z)] Jo ((1 - z)bq) fs{x, 6, q) 



and similarly in the integrand of ()67|) the term 

a,(g2) 



[2;P,,(2;) + zPggiz)] Jo ((1 - fcix, b, q) 



Thus we obtain 



X <^ Jo((l-^)g6) 



g2 27r 

eiz - x) [P,g{z) fs b, g) + P,,(^) fc 6, g 



- [zPgq{z) + zPgg{z)] fs{x, b, q) 



- [zPgg{z) + ZPgg{z)] [l " Jq ( ( 1 - Z)?^)] ^(X, 6, ?) 



(68) 
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X !.M{l-z)qb) 



Q{z - x) P,,{z) fs ( ^, b, g) + Pgg{z) fc b, q 



- l^Pggiz) + zPqg{z)] /g(x, b, q) 



[zPgg{z) + zP,g{z)] [l - Jq ( ( 1 " /^(x, 6, g) 



(69) 



We consider this equation in the region 6 ^ l/Qo^ set 6 = in the inhomogeneous term, 
make the approximation {1 — Jo((l ^ z)^!^))} — — z)qb — 1) in the last terms and hmit 
the integration to g < ■min{Q, 1/6} in the remaining terms, to get for 6 > 1/Q : 

fs{x,b,Q) = ^xqsix,Ql) + I j\z'^^Q{q'~ql)i^Q{l/b-q) 
Q{z - x) [P,,{z) fs (^, 6, q) + P,g{z) fa (^, 6, q 



[zPqq{z) + zPgq{z)] fs{x, 6, q) 



e(g - q) Q(q - 1/b) 6(1 - l/(g6) - z) [zPqq{z) + zPgq{z)]fs{x, b, q) (70) 



fG{x,b,Q) = -xg{x,Ql 



X 



9(2; - x) Pgg{z) fs ( |, b, g) + P33(z) /g 6, g 



+ ^Pqgi^)] fcix, b, q) 



e(Q _ q) Q^q _ 1/5) e(l - l/(g6) - z) [zPgg{z) + zPqg{z)] /g(x, 6, g) \ (71) 



where 



qs{x,Ql) = ^{qi{x,Ql) +qi{x,Ql)) . 



(72) 



i=l 



We note that fi{x, b, q) and fi{x/z, b, q) in all terms, except those in the last terms of eqs. (fTOj) 
and (f7T|) . are in the region g < 1/b. We shall show that those terms give l/2xqs{x, 1/6^) and 
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l/2xg{x, l/&^), respectively. Let us observe at first that for Q < 1/6 eqs. (fTOI) and (f7T|) read: 



- [zPqg{z) + zPgq{z)] fs{x, b, q) 



X \e{z-x) P,,{z)fs[^,b,q)+Pgg{z)fG[^,b,q 



(73) 



- [zPgg{z) + ZPqg{z)] f g{,X , 6, ^ , 



(74) 



wliicli are just tlie DGLAP equations in tlie integral form for l/2xg{x, Q"^) and 1/2 xqs{x^ Q"^). 
We therefore get fc^x^b) = l/2xg{x,Q'^) for b < 1/Q. The first terms in eq. (fTUj) and (f7T|) are 
the same as in the right hand side of equations (fTHj) and (fTIj) except that instead of Q{Q — q) 
we have 9(1/6 — q). The first terms in eqs. (fTOj) and (jTH) reduce to 1/2 xg5(x, 1/6^) and 
l/2xg{x,l/b'^). Therefore, eqs. (ffn|) and (ffT|) read: 



fs{x,b,Q) = -xqs{x,l/b^) - / dz 



dq^ as{q^) 
g2 27r 



X e(g - 1/6) 0(g - q) 6(1 - l/(g6) - z) [zP,,{z) + zPg,{z)] fsix, 6, q) (75) 

dq"^ a,(g2) 



fG{x,b,Q) = ^xg{x,l/b'^) - j dz j 



g2 2ti 



X e(g-l/6)e(Q-g)e(l-l/(g6)-^)[zP,,(z) + zP,,(z)]/G(x,6,g), (76) 
and their solution read: 



fs{xAQ) = -exp{l/2 5,(6,Q)}a;g5(x,l/62) 



/g(x,6,Q) = -exp{l/2 5G(6,g)}a;^7(a;,l/62) 



with Sq and defined by: 



S,{b, Q) = -2 I 4^ / dz [zP,,{z) + zPgqiz)] 



SGib,Q) 



i/fe2 27r Jo 
"^"dq^asiq") '^'^'^''^ 



1/62 27r JO 
15 



dz [zPgg{z) + ^Pgg(^)] 



(77) 
(78) 

(79) 
(80) 



Substituting explicit expressions for the splitting functions Pij{z) and setting 1 instead of 
1 — 1/ (bq) as the upper integration limit in the integrals over non-singular at z = 1 terms, and 
using eqs. (fT ^ - (fT^ we get: 

Saib.Q) - - r WS H^b^) + 4"1 (81) 

with A^Q and Bq^ given by eq. (fTTj) . We also find that Sq{h, Q) is given by eq. As in the 
case of non-singlet quarks the first term on the right hand side of eq. ()81|) . with the integrand 
proportional to ln(g^6^), gives the result which is given by the term with Aciots) in eq. (jH)) and 
the fixed coupling only. Thus, the subleading differences are due to a different scale of the QCD 
coupling. We would get exactly the same result for this term choosing the scale of equal to 
[g(l — zyf' instead of q^ in the integrand defining l\{h,Q), eq. (pT|) . The second term on the 
right hand side of eq. (jHTj) is identical to that defined by Bg{ols) in eq. As in the case of the 
valence quark distributions also for the gluon and sea quark distributions the non-perturbative 
effects will enter as the Q— independent factor(s). 



6 Numerical results 

We have shown in the previous sections that the parton shower described by the CCFM equa- 
tions in the single loop approximation (PT|) - (jin|) generate complete resummation effects in LO, 
i.e. it automatically resums the double and single "large" logarithms \yi'{Q'^IP') in the region 
^ 1. In this Section we present numerical results based on the exact (numerical) solution 
of the CCFM equations in the 6— representation and compare them with the predictions based 
on the resummation formulas. In our analysis we assume a factorisable form of the initial 
conditions (jl^ 

= ]^F{h)xpu{x,Ql), (82) 

assuming for simplicity the same input profile F{h) for quarks and gluons. It is determined by 
a non-perturbative fcr— distribution at the scale Qo, and at 6 = we have the normalisation 
condition F(0) = 1. Thus, we assume the Gaussian form of the profile function 

F{h) = exp{-6V6^| ^ (83) 

where we set 6q = 4GeV~^. The starting integrated distributions pk{x,Ql) are taken from 
the LO GRV analysis 36 . With the factorisable form (jH^ of the starting quark and gluon 
distributions, we obtain the corresponding factorisation of the function (j^Uj) 

Wij{xi,X2,b,Q) = F^{b)Wij{xi,X2,b,Q) , i,j = q,q,g. (84) 

For the soft gluon resummation, the expression for the function Wij takes the following form 

W~{x,,X2,b,Q) = exp{l/2[S,{b,Q)+S,ib,Q)]}x,p,ix,,l/b') X2P,(x2, 1/&') . (85) 
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In Fig. 1 we show results for Wgg{xi, X2, b,Q = M^) plotted as the function of 6^, calculated 
for xi = X2 = Mh/^ for equal to 1.8 TeV (Fig. la) and 14 TeV (Fig. lb). The Higgs 
boson mass was assumed to be equal to 115 GeV. We show the results based upon the numerical 
solution of the CCFM equations (PT|) -(pS j) (sohd hues) and confront them with the predictions 
based upon the soft gluon resummation formula with the argument of the exponent 5'g(6, Q) 
given by eq. (jHT|) (dotted lines). We also show the results corresponding to Scip, Q) given by the 
standard eq. © (dashed lines). In both latter cases the scale was replaced by Q"^ / (fo^Q^ + l) 
that guarantees smooth transition to the region IP' < 1/Q^. We have also replaced Scip, Q) by 
S^J^{b,Q), setting S^J^{h,Q) = in the region where 5^(6, Q) > 0, and S^J^{b,Q) = Saib^Q) 
in the region where S'g(&, Q) < 0. We can see that the function Wgg calculated from the 
exact solution of the CCFM equation stays most of the time below the predictions based upon 
the resummation formulas. Only at relatively large values of 6 ~ 1 GeV~^, the resummation 
formula corresponding to S'g(&, Q) given by eq. (|HT|) leads to smaller magnitude of Wgg than 
that corresponding to the exact solution of the CCFM equation. This difference is caused by a 
different behaviour of the Sudakov-like form factors Tg = exp{SG{b,Q)} in all three cases, as 
can be seen in Fig. 2. In this figure we plot as the solid line the function exp{S'g'"^^^(6, Q)} 
with 

S%^^^\b,Q)= / dzz[Pgg{z) + P,g{z)]{U{l-z)qb)-l}, (86) 

Jql T ^TT Jo 

and compare it with the functions resulting from the soft resummation functions with Sq given 
by eq. (jHT|) (dotted line) and eq. (jHI) (dashed line). 

In Fig. 3 we show the impact of those differences on the transverse momentum distribution 
of Higgs bosons produced at ?/ = 0, eq. (P). We present the results for ^/s = 1.8 TeV and 
for y/s = 14 TeV. We can see that the results based upon the exact solution of the CCFM 
equation (solid lines) gives smaller magnitude of the cross section at its maximum that those 
generated by the soft gluon resummation formulas (jHT|l (dotted lines) and (dashed lines). 
The former tends, however, to generate the longer tail of the pt distribution. 

7 Summary and conclusions 

In this paper we have discussed the transverse momentum distributions of the Higgs bosons 
within the unintegrated distributions framework. We considered the gg ^ H fusion process 
which is the dominant mechanism of inclusive Higgs boson production. The unintegrated 
distributions were obtained from the CCFM equations in the single loop approximation which 
is equivalent to the LO DGLAP evolution for the integrated distributions. We have utilised 
the fact that the CCFM equations in the single loop approximation can be diagonalised in the 
b representation. We have shown that the conventional soft gluon resummation formulas are 
embodied in the solution(s) of the CCFM equations. We have examined possible differences 
between the exact (numerical) solution of the CCFM equation and their approximate solution 
leading to conventional soft gluon resummation formulas. 
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Figure 1: The function WGG{xi,X2,b,Q) as a function of b for the Higgs boson production 
at Q = rriH = 115 GeV in the gluon - gluon fusion at the Tevatron (a) and the LHC energy 
(b). The sohd hues correspond to the resuh with the unintegrated gluon distributions from the 
CCFM equations. The results with the LO soft gluon resumation formulas are shown by the 
dotted lines {So given by eq. (jHH)), and the dashed lines (5*^ given by eq. Q). 



21 



,0 

10^ 10" lO'' id" 10' 10° io' 

b' [GeV "'] 



Figure 2: The Sudakow-like form factor Tq = explSc}, as a function of the impact parameter 
b, from the CCFM equations (sohd hne) and the soft gluon resummation formulas (dotted and 
dashed hues). The function Sq is given by eq. (jHBj) for the sohd hne, eq. (jHT|) for the dotted 
hne, and eq. Q for the dashed hne. 
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Figure 3: The transverse momentum distribution of Higgs boson from the gluon-gluon fusion 
at Q = niH = 115 GeV for central rapidity, y = 0. Figures (a) and (b) correspond to the 
Tevatron and LHC energy, respectively. The solid lines are based on the numerical solution of 
the CCFM equations PTj) - (P5jl . while the the dotted and dashed lines results from the LO soft 
gluon resummation formulas with So given by eq. (|HT|) and eq. 0, respectively. 
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